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CN ■ Abstract 

' We show that, in quadratic lagrangian theories of gravity, isotropic cosmological singularities are stable 

^ I to the presence of small scalar, vector and tensor inhomogeneities. Unlike in general relativity, a particular 

exact isotropic solution is shown to be the stable attractor on approach to the initial cosmological 
' singularity. This solution is also known to act as an attractor in Bianchi universes of types I, II and 

IX, and the results of this paper reinforce the hypothesis that small inhomogeneous and anisotropic 
perturbations of this attractor form part of the general cosmological solution to the field equations of 
. quadratic gravity. Implications for the existence of a 'gravitational entropy' are also discussed. 

00 . 

^ ■ 1 Introduction 

(N 
O. 

, There has been considerable interest in the stability properties of isotropic initial singularities in general 

' relativistic cosmologies. At first, this question was bound up with the issue of whether the initial singularity 

, itself was a robust prediction of general relativistic cosmologies in the presence of typical fluids, like dust and 

radiation, or whether it was merely an aberration of poor coordinate choice or high symmetry. This question 
^ , was answered unambiguously by the proofs of the first singularity theorems which approached the problem 

by means of the geodesic equations and the causal structure of space-time instead of the Einstein equations 
[1]. However, there still remained the question of whether any initial singularity was likely to be isotropic. 
The discovery in 1967 of the high level of temperature anisotropy in the microwave background made this 
question one of central importance for observational cosmology [2] . It turned out that the general behaviour 
' of cosmological singularities dominated by matter or radiation was highly anisotropic [3, 4], even chaotic 

[5] , and any attempt to explain the existence and rotation of galaxies by primordial turbulence required the 
initial singularity to be highly anisotropic [6], while attempts to dissipate strong initial anisotropics in the 
spirit of the chaotic cosmology programme [7] had to face the problem of entropy over-production in the 
very early universe [8]. However, it was possible for singularities dominated by stiff matter to be close to 
isotropy and 'quiescent' because the energy density of this extreme form of matter, characteristic of a massless 
scalar field, diverges as fast as the most rapidly diverging anisotropy modes on approach to a cosmological 
singularity [9, 10]. More recently, the advent of inflation [11] as a potential explanation for the observed 
isotropy of the visible part of the universe removed the need to appeal to a globally isotropic singularity in 
order to explain astronomical observations, or for anisotropy to be dissipated [12]. It is sufficient to explain 
the existence of local isotropy on the scale of the observable universe and this is achieved by natural means 
in the chaotic inflationary scenario [13]. Despite the success of this approach, there remain a number of 
questions whose answers depend upon the stability properties of an initial isotropic cosmological singularity. 
A sufficiently isotropic and homogeneous horizon-sized patch is needed somewhere in order for infiation to 
occur, and attempts to understand the quantum nature of gravitation will need to understand whether a 
physical cosmological singularity generically occurs and, if so, what its most likely properties are. 
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Penrose [14] has proposed that initial c;osniological singularities are not simply signals that the underlying 
gravity theory has broken down, but play a crucial role in defining the arrow of time. Initial singularities 
are thereby required to be highly isotropic - so possessing low gravitational 'entropy', perhaps identified 
with the Wcyl curvature - but any final singularity in a collapsing closed universe would necessarily be 
strongly anisotropic, with high 'gravitational entropy'. Since such isotropic singularities are special in general 
relativity [15], it is necessary to introduce this gravitational entropy hypothesis as a boundary condition in 
order to justify choosing a highly isotropic initial singularity. 

Other boundary conditions can also be chosen which make the present level of expansion isotropy un- 
derstandable even if the initial state is not isotropic. Barrow [16] showed that the natural requirement that 
all energy densities (including that in anisotropising gravitational waves) are smaller than, or of order, the 
Planck density at the Planck time leads to a general requirement that the microwave background anisotropy 
be no larger than 0(10"'') today. Typically, inflation scenarios assume that arbitrarily large energy densities 
do not arise [13]. 

Here, we want to explore the stability of isotropic singularities in higher-order generalisations of general 
relativity, in which the Einstein-Hilbert lagrangian is supplemented by the addition of the quadratic scalar, 
R, and Ricci, RatR"''', curvature invariants. Since initial cosmological singularities arc expected to involve 
infinities in one of more of the curvature invariants of the space-time, we expect that the addition of higher- 
order terms to the lagrangian would produce a new dominant behaviour to such singularities. We would 
expect the higher-order terms to control the dynamics as the singularity is approached, and to determine 
the stability properties of any isotropic special solution in that limit. In this paper we are going to consider 
the contributions of the quadratic Ricci terms, RabR"^, to the lagrangian. The quadratic scalar, R^, contri- 
butions are conformally equivalent to the presence of a self-interacting scalar field and are understood. In a 
subsequent paper, we will study how our results may be modified by complicated effects that arise at general 
order, when (RabR"^)" terms are added to the lagrangian. 

Specifically, we will investigate the stability of an initial isotropic zero-curvature Fricdmann singularity 
to small inhomogeneous scalar, vector, and tensor perturbations in generalisations of general relativity which 
contain quadratic Ricci and scalar curvature contributions to the gravitational lagrangian of general relativity. 
We will show that there exists a particular isotropic solution of the pure quadratic Ricci theory which is a 
stable attractor as t ^ for small inhomogeneous perturbations of the metric when the gravitation theory 
created by the addition of quadratic Ricci terms to the general relativity lagrangian. Previous studies have 
shown this particular isotropic solution to be a stable f, ^ attractor for spatially homogeneous vacuum 
universes of Bianchi types /, //, [17] and IX [18]. Our results lend further weight to the conjecture of Barrow 
and Hervik [17] that this isotropic vacuum solution of the pure quadratic Ricci theory characterises part of 
the general cosmological solution in these particular higher-order generalisations of general relativity. This 
behaviour is completely different to that found in general relativity, where isotropy is unstable in the t ^ 
limit of an initial cosmological singularity [3, 4]. 



2 A Special Solution 



Consider a quadratic gravity theory with action 




The field equations obtained by varying this action with respect to the metric are [20, 17, 21]: 



'6> 



a 



(1) 



where 



p(2)a _ _1 

' ~ 2 



R^g^ + 2RRI + 2glUR - 2i?;"b , 
R'"^Rcdgb + -I- -g^HR — — 2R°^cdbR' 
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where Gab = Rab — \Rgab IS the usual Einstein tensor, is the energy-momentum tensor of the matter 
sources, and A is the cosmological constant^. 

We take A = and consider perturbations about the spatially flat, homogeneous and isotropic FRW 
spacetime, 

ds^ = -dt^ + a^(t){dx^ + dy^ + dz^), (2) 

with scale factor a(i)and associated Hubble expansion rate H = d/a, for which the background equations 
are: 

i/K = + 2{B + 3A) (2HH - + 6H^H^ , (3) 

-lif^+P) = H + 2{B + 3A)(h + 3HH + 6H^^, (4) 

where we have defined a density, fi = —Tq, and pressure p= ^T^- 

For radiation obeying the equation of state p = |/U, there exists a special exact solution of these equations 
with [19, 21] 

a{t) = ti. 

Importantly, this is also an exact vacuum solution (/U = p = 0) of the purely quadratic theory^ {B + 3A^ 
oo), i.e. it solves 

2HH -H^ + 6H'^H = 0. (5) 

Integrating this equation, we have 

with M constant. Since H = we can take the square root and integrate to obtain 

3 3a3' 

with D constant. Hence, we can see that a ^ as t ^ and const +0{a~^) to leading order as 

a ^ oo. So the purely quadratic vacuum solution mimics radiation early on, then dust plus a cosmological 
constant at late times. 

Clifton and Barrow [23] found that for a universe filled with a fluid with equation of state p = uj^i in a 
theory where the lagrangian depends only on the Ricci-squared term like {RabR""^)^ , there is an isotropic 
fluid-filled FRW solution with 

a{t) = t^i^, 

and also that there exists an isotropic vacuum solution with 



a = tV2 :n=l, (6) 
a = :ny^l, (7) 



where [23] 



3(1 - 3n + 4^2) ± ^^3(48714 - AOn^ - bn^ + lOn - 1) 
" 6(1 -n) ■ 



^For a different covariant formulation see ref. [22]. 

^Barrow and Hervik [17] showed that this is also true for the case of Priedmann models with non-zero curvature, in which 
case the curved Priedmann radiation solution is an exact solution of the full equations in the presence of radiation and also of 
the purely quadratic vacuum equations (B + 3j4 — > oo). 
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In particular, for the case of n = 1 and w < 1, we see that 



> I > |, so the special vacuum solution 



3(w + l) 

will dominate over the fluid- filled effects in the field equations ast^O. 

The special isotropic vacuum solution (6) is of special interest. Previous studies have shown that it is 
an attractor for solutions of Bianchi types I and II as t — > 0, [17]. It is also a stable attractor in Bianchi 
type IX 'Mixmaster' universes as t ^ [18]. This is in complete contrast to the situation that exists in 
general relativity (A = B = 0) where the shear and curvature anisotropy terms dominate the dynamics as 
the singularity is approached at t = 0, [3, 4]: in the case of Bianchi type IX, that dynamical approach is 
chaotic [5]. As a result an isotropic initial singularity constitutes a set of zero measure in the initial data 
space. In quadratic gravity, the situation is completely different, the results in Bianchi types I, II, and IX 
suggest that the special isotropic solution may be an attractor for the time evolution of a part of the general 
solution of the cosmological evolution equations. In the sections to follow we will show that this special 
vacuum solution is also stable against the effects of small inhomogeneous perturbations (scalar, vector, and 
tensor) as t ^ 0. 



3 Inhomogeneous Perturbations 

The general perturbation of the FRW metric may be written as 

ds^ = -a^(l + 2a)dr]^ - a^B^drjdx" + 0^(6^/3 + Cai3)dx"dx'^ , 



(8) 



where is a conformal time coordinate that is related to the comoving proper time, t, by dt = adr/. We can 
decompose the perturbation variables into their "scalar", "vector" and "tensor" parts, as follows^: 

Ca/3 = 20(5^/3 + 27 + 2C(a,l3) + 2Ca/3. 

There are four scalar perturbation variables, a,/3,<p and 7, two vector variables, B^ and C„, and one 
tensor perturbation, Cap. The quantities B^ and Cq, are divergence- free, i.e. B",a = = C"q,, and Cap 
is transverse and tracefree, i.e. C"/3,a = = C"c«- These three types of perturbation evolve independently 
of each other at linear order. The linearised equations governing the perturbations in the general quadratic 
gravity are presented in a scries of papers by Noh and Hwang [24, 25, 26]. 

For the gravitational-wave (tensor) perturbation modes in the general quadratic gravity, the equation for 
the perturbation is ([24]): 



STS 



D'^ + 2A (rDI + RC^) - + 3HD'^ _ ^ (z)^ + 4HC^ 

(ii + H^^ D'^+(h + HHj Cf, 



where 



For the vortical (vector) perturbations, following Noh and Hwang [25], we set Ba{t,x) = b{t)Ya{yi) , with 
AYa = —k^Ya and use the gauge-invariant variable '^Ya = B^ + aCa- In the gauge where = 0, the 
equations are 



2 k"^ 
* il* - I -i? - ^ 1 * 



-,^il + 2AR)^-B 
and the conservation of angular momentum transparently holds: 

a'^{fj,+p)vu = O 



= {lJ'+P)Vu. 



^Note that there is a factor of 2 difference in the d'qdx" component from the metrics given by Noh and Hwang in [25] and 
[26] . Defining poa = ~<J^ {/9 « + ) consitent with their equations and appendix quantities. 
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with Q constant. 

Finally, for the scalar density perturbations, we introduce the combinations x = + 
Z{Ha — (j)) — Then we have the following set of equations ([26]): 
Energy constraint: 



(5r° = 2 



A(j)\- 



A<j> 



Hk] +R{k + ma) - iHSR 



3{H + 



A 



6R 



a? J 



+ 



4:{H - QH'^) 



2A 



Hk 



A(j) 



+ 3HdR 



- 3H^ + 



SR + 6HHa + 6 -HH + 2H^ - 9H^H a 



8H A 



6H + 5H'] n-^=r 



Momentum constraint: 



3a 
-B 



Ax 



Ax 



3 a2 



A -2Rk - 2R 



Ax 



3Ra + 3SR - 3HSR 



3H 



A 



12H' + _ K + ^ 



Ax\ 3SR 3H5R 



2 +■ 2 



-3Ha +^(2H + 9H'^y a + ^ 2iJ^ + (i? - 2H'^^ 
Tracefree propagation: 

5T^-^dT^6^ = 1 (^V«V^ - i<5^A^ x + Hx- a + 'iA(^Rx+ {HR + R)x- SR- ^- a) 



B\{x + Hx-ct>- a)- -H{x + Hx-ct>- a)' 



8(£r + H^) + — 



{X + Hx - (i>- a) 



+5R-m4>-&{H + H^)x + ^-[^ + ^)] 



Raychaudhuri equation: 



5T:t-5T^ = 2k + mK + 2(3H + ^a + 2Al^Rk+(AHR+R^K 

A 



6i? + 3HR + 2R[3H + 



AHR+R] K + 3Ra 



a-3SR- 3H5R +{6H^ + 



A 



SR 



-B{ -A[Hk + 



+ 



A(j)\- 



12H Hk 



A(j)\ 



-8(h-6H^)+^'\ (HK+^y25R 
+6H'^SR - 28Hk - 6 (2H + 5H'^^ a -12 (2H + 5H'^^ " + 3H'^H 



a 



16H 
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Trace equation: 



and 



5T = -SR-2{SA + B) 



5R 



5R + ■^H5R - ^6R - R{k + a) - (2R + 3HRj a 



A 



Note that no choice of gauge has yet been made for the scalar perturbations; we will ultimately want to 
work in the gauge 5R = 0, so that wc arc considering perturbations in flat space. 

We are primarily interested in the effects of the RabR"''' term and so henceforth we set ^ = 0"^. We 
consider the three perturbation modes in turn in order to determine whether the isotropic and homogeneous 
vacuum solution (6) is stable against their effects as t — > 0. 



4 Gravitational-wave perturbations 

The metric for the gravitational-wave perturbations takes the form 

ds^ = -dt^ + a^{Sai3 + 2Cocf))dx"dx'^, 

C^ = ^ = Cl^ (9) 

where (7^/3 = Cai3{x,t) and is gauge- invariant. For convenience, from now on we drop the indices on C^/j . 
The remaining perturbation equations in vacuum are: 



= D-B{D + SHD-^{D + AHC)-6 {H + H^)D + {H + H^)C 



D = C + 3HC-^, 

where £» = D| and C = C|. 

Substituting H = l/2t, and taking the B — > oo limit, keeping only {} = 0, this reduces to 







\ 




- 








T 




AC 


D. 






2t 







(11) 

Substituting for D using (11), we have the single 4th-order evolution equation for the metric perturbation 

C: 

„ .... 3... 3 ^ 2AC AC AAC 

4.1 Large scales 

In the long-wavelength limit, the last three terms in (12) can be neglected and we have 

o=c + lc + ^c, 



*Note that for our special solution, eq. (6) with a{t) = , wc have i? = to background order. Consequently, the perturbed 
parts of the gravitational wave and vorticity equations arc independent of A. This is also true for the scalar perturbations when 
we choose the uniform curvature gauge, i.e. SR = 0. So our results are also valid for the general quadratic gravity. 
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and so the solution is a linear sum of power-laws with 



hence 



C oc 

= g(g-l)(2g-l)(2g-3), 



C = £ + At + + (t2 ^ 



where e, A, v, are independent arbitrary tensor functions of x satisfying the same constraints as C, given 

by (9). 

Note that as t ^ 0, the perturbation C const and the a — solution is stable against gravitational- 
wave modes. The leading-order terms in the metric perturbation expansion are therefore (restoring the tensor 
indices) : 

ds'^ = -dt^ + t{5o,i3 + 2ec«/3 + Si^a/?*^ + "^Kpt + 2C ^ fjti ) dx" dx/^ , 

where the ^^/^(x), z/a/j(x), Aq/3(x), ^^^(x) are arbitrary traceless symmetric 3x3 tensors subject only to (9) 
and coordinate transformations and the constraint equations. 

4.2 Small scales 

If we look for a separable solution to (12) of the form 

C = T{t)X{x) 

then 

.... 3... 3 ^ 2AC AC AAC 

C+tC + TT^C' = —— + -5 7T- 

3... 3 •• 2fAX TAX TAAX 
=> TX+-TX + -.TX = ^— + —2 

If AX = ex, where 6 is constant so A AX = 6 AX = e'^X, then 

.... 3... 3 " 2fe 1 6 e^T 

f f + St f + (^^ - 2^^*^ f - f + 9^T = 0. (14) 

li 9 = 0, then we obtain the same results as in large-scale limit. If we take the late-time approximation 
or the small-scale limit by keeping only the terms in < in (13) then we need to solve 

2tf + f-0T = 0, 

so as 6 will be negative, 

T = Ci cos{V-2et} + Ca sin{v'-26lt} 
which approaches a constant as t ^ 0. 
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4.3 The general sepeirable solution 

If we change to conformal time, ^ = X' P^* ^ ~ then the master equation (14) is transformed to 



T"" + 2fc^T" + fc*T = 0. 



The general solution is 



T*{t) = Ai cos[2kVi] + A2 sin[2A;v^] + AaVi cos[2kVi] + A4Vi sm[2kVi] k ^ 0, 

T = Bi+ BiVt + Sgt + BaIVI k = 0. 



We note that, for A; ^ 0, 



T*{t) ^ Ai + {2kA2 + A3)Vi + 0{t) const 



as t —^ 0. 

Hence the full separable solution for the metric perturbation is 

AX^0 = -fc2x„ 



5 Vortical Perturbations 

The metric expansion for the vorticity perturbations takes the form 

ds"^ = -a^drf' - 2a^ Bocdrjdx"' + {5ai3 + 2C(a,/3)) dx°'dx^, 



(15) 



with = = C«„. 

Following Noh and Hwang [25], we set Ba{t,x) = b{t)Ya{x) , with AY^ = —k'^Ya, and use the gauge- 
invariant variable: 

^Yoc = Ba + aCoc. 

In the gauge defined by Cq, = 0, linearising about the special vacuum solution with H = ^^and R = 
we then have 



2a2 



* + — H * 

2t t 



and the conservation of angular momentum gives: 
with O constant. 

5.1 Lcirge scales 

On large scales we have 



2t 



2n 



In the case where the (quadratic) B term dominates on the LHS, we have 



(16) 



and ^ — > const as i ^ and the isotropic metric is stable. 
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If we keep the general-relativistic ^ term on the LHS of (16), and write ^ = tif, then we have 



1 \ „ 20 



3 



This is an inhomogeneous modified Bessel equation; so, if we set 

it 1 
y= ^'^= 7^9 



we obtain 

y'f" + yf + {y'-'^')f = 9ix) 

and its solution is 



-IjJ^) f^Yj^] MS^]^dt. 



2 3 V/By J y/B ^ \y/Bj \ Bk"^ J ^ 

Hence, on large-scales the general solution in the Cq, = gauge is 

= *F„=t3/(f)F„(x), (18) 
AYc = -k^Yoc. (19) 

Again we note that B^ — > constant as t — > and the isotropic solution (6) is stable in that hmit. 
5.2 The pure quadratic theory 

Now consider the case where the B-term dominates, but this time without taking the long-wavelength limit, 
i.e. we drop just the GR term on the left hand side. In that case the governing equation is 

''+2t^—=-BkH 

and so 

2f] 

* = £>i cosh{2kVi} + £>2 sinh{2fc\/i} - (21) 

where Di and D2 are constants. 

Since B^ = ^'^^(x) in our chosen gauge, we again see that B^ — > constant as t — > and (6) is stable in 
that limit. 

6 Scalar Perturbations 

For the scalar perturbations, the metric will take the form 

rfs^ = -a^(l + 2a)dr)'^ - 2a'^ l3 ^^drjdx"' + a^(5a/3(l + 2^!>) + 27_„^)da;"rfa;^. 
There are possible gauge choices that can still be made and we shall consider two useful choices. 
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6.1 Conformal Newtonian Gauge (x = 0) 

The gauge choice x = o,(P + o-i) = (zero shear) greatly simpUfies the perturbation equations, since we just 
have two perturbed variables to worry about. They reduce to the following set: 
Energy constraint: 



ST^ = -B<i 6^ + 15^ -6^ -4^ -12^ + 2^ + 2 



+ 



2A(/) 



3(/) 

T 



i2 

3a 
2*2' 



T 



a 3d Aa A^c 

+ 2^3 +l3-+2^ 



(22) 



Momentum constraint: 



T° = -—V 

-'a ,1 
t2 



B 



„ „ fli A(/) „ Afl!) a , fi Ad Aa 



a 
2t 



(23) 



Trace-free propagation: 



5T^--ST:^5% = - 



i(v«V,-i.,"A) 



B<b(b+ — 3 OL — 

^ ^ 2t t 2t t 



+ + a 



(24) 



Trace equation: 



5T 



= -2b|6'^" 



6 



2Aa 



10 



A^ A<j> Acj) 



t 



10 



■■ 12<P 4A0 3d 

—od) h 

^ t t t t 



,A2 



3a 3d Aa 2A2a 



(25) 



6.1.1 Large-scale limit of the pure quadratic theory {B oo) in the zero-shear gauge 

If we ignore the effects of matter^, i.e. 6Tl^ = 0, and take the large-scale limit then each term in curly 
brackets vanishes. The combination j (23) + | (22) — ^ (24) gives 



2^ + 6^ + ^ 



(26) 



to which the solution is 



Using the energy equation, we then obtain 



a = Cit-^ + 3(72*5 + Csi'. 

So, at first sight it appears that there may be an instability as t ^ 0; however, we note that the coefficients 
of the < 2 terms in a and tj) are the same. This means that they can both be locally defined away by the freedom 
to choose our initial time. If we set t = f — fo, then a(t) = 1 5 = (t — to)'^ = i^{l — toi~^)^ i's — ^toi~2 . In 
fact, by choosing a different gauge, we see that this is just a curvature perturbation when we compare our 
flat FRW background with the similar open and closed models. 



oc a-3(i+") oc t- 



for radiation, whereas on the right-hand side of the energy equation typical terms are ' 
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6.2 Uniform Curvature Gauge (SR = 0) 

If wc use the uniform curvature gauge then i? = and we are perturbing in flat space. The set of perturbations 
equations are now: 
Energy constraint: 



= -3- + 2 



A2 



A(p 3a Ax 



t 



+ 777^ - 



2t2 t2 



3(1) 

2*2 



A0 A(f) 



+ 21-V + 2— +6^-7 



"^3 



t 



t2 



2-^ + 777^7 



3a 15d 3a 7Aa Ax 3Ax HAx A^x 



t2 



2*2 4*3 2*4 2*3 



*2 + 2*3 



+ 



2*4 



*3 



(27) 



Momentum constraint: 



Q = 3<i)- 



3a 
27 



rJv;^-^^'^ q*^ A0 A0 3d 9d 3Aa Ax 



(28) 



Trace-free propagation: 







7(j) A(j) a Aa 



3x 3x Ax Ax 
* 2*2 



^ + 4*2 4*3 



^-TT^ . (29) 



The trace equation is trivial, but we have a fourth equation from 5R = 0: 

= -3^ - 6 



^ A^ 3d Aa_ Ax_ Ax 
* * 2* * * *2 ■ 



(30) 



We can now eliminate x from our equations (from the momentum constraint, we see that x is stable 
whenever a and (j) are). 

Using (28), equation (30) becomes: 



= 



B 



+ 4*20 - 



*2d 3*a 



- *' 



.3' J,' 
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(31) 



+*^A^ + 2*A(/> - 



*2d 3*d *Ad 



Equation (27) becomes: 
= -b\2 



A0 A</. A2^ 
2* + *2 + * 



A2 



+ 



-*( 



36 ^6 Ad) a 3d 
_z: + 2--^ — --i 

2 * 



* 



4* 



a 1 
2*2 B 



*0- 



*2 *2 2*3 2*4 

a 



~2t 4*2 



a 



2*2 



(32) 



6.2.1 Large-scale limit of the pure quadratic theory {B — > oo) in the uniform curvature gauge 

In the large-scale limit of (31), setting (j) = E\t^, a = Fxt^ and dropping the B~^ terms, we find that 



1 



A(A-l) A-- (A+l)i;A--^ 



Fx 



0, 

and in the large-scale limit of (32), keeping only the terms inside {}, we have 

3' 



(A - 1) £;aA(2A - 1) - Fa a 



0. 



(33) 



(34) 
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The right-hand bracket of each of these last two equations cannot simultaneously vanish unless E\ = 
Fa = 0, so A = 0, or 1 and Fq = Fi = 0. There is a further restriction on Ei and Fi from (29), specifically 



Hence, we have the solution 



(j) = Eo + Eit^ +Eit, (35) 

g 

a = -Eit, (36) 



\x = ^{t^fi + ti) = -lE.t^-^E^t. (37) 

Assuming Ax = —k?x, we have 

/3 = 2^^i-*^5(i,x), 

There is a remaining freedom in our gauge choice. In the conformal Newtonian gauge, it is usually used 
to set both /3 = and 7 = 0. The other parameters are unaffected by this choice. Here, we can use this 

freedom to set g{t,x) = 0. 

Again, we see that the metric perturbations approach constant values as i ^ and the isotropic solution 
is stable against scalar perturbations in this limit. 

7 Summary 

Combining our results from the quadratic theory, to leading order in the time as f — > 0, we have found that 
the metric in the neighbourhood of the special solution (6) with isotropic scale factor evolution a = t^^^ 
takes the general form 

ds'^ = -0^(1 + 2a)dr]^ - a^Badrjdx" + a^iS^/} + Cap)dx'^dx'^ , (38) 

where 



a = -Eit, (39) 
5 

B„ = ■^Ei^^ + 2lDi-^ + 2kD2t-^+2k^DitjYa{x), (40) 

Co,f3 = 2(^Eo + E,t^ +Eitya0 + ^Ei,c,fs+ (41) 
+2 (Ai + {2jA2 + A3)t^ + 2j{jAi + A4)t) X„^(x). (42) 

where AYq, = —k^Ya, AXap = —j'^Xap, the Ai and Di are constants, and Ei = Ei{x). 

These results are striking. The special isotropic and homogeneous solution with scale factor evolution 
a — t^/"^ is found to be stable against the effects of small scalar, vector and tensor perturbations in the t ^ {) 
limit. This adds to the evidence gathered from earlier studies of the stability of this solution in the Bianchi 
type I, II and IX universes [17, 18]. It adds further weight to the conjecture that small perturbations of 
this isotropic singularity form part of the general cosmological solution for theories of gravity with quadratic 
lagrangians of the sort studied here. This behaviour is quite different to that found in cosmological solutions 
of general relativity, where isotropic expansion is unstable as t ^ Q. It arises because the higher-order 
quadratic terms in the lagrangian contribute terms which dominate on approach to an initial singularity 
(typically as 0{t~'^) compared to general relativistic contributions at 0{t~'^)) contribute isotropising stresses 



12 



which dominate over the shear stresses arising from expansion and 3-curvature anisotropy, which dominate 
the initial singularity in general relativistic cosmologies. 

These results have a number of implications for proposals to introduce special cosmological initial condi- 
tions in order to ensure that the initial state is isotropic, with low gravitational entropy. It appears to make 
additional stipulations of special initial conditions unnecessary in quadratic gravity. However, it may cause 
problems for the gravitational entropy scenario in closed recoUapsing universes because we expect the same 
quadratic stresses also to drive the solution towards isotropy on approach to any future 'big crunch' singu- 
larity. This would prevent it being the high gravitational entropy state that a 'Second Law' of gravitational 
thermodynamics would lead us to expect. 

Finally, despite the special features of the cosmological models in quadratic theories of gravity, we need 
to understand if these unusual results regarding the stability of isotropy persist when higher-than-quadratic- 
order contributions to the gravitational lagrangian are included. In a subsequent paper we will present the 
results of this more complicated study. 
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